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Abstract 

Examples of symmetric informationally complete positive operator valued mea- 
sures (SIC-POVMs) have been constructed in every dimension < 67. However, 
it remains an open question whether they exist in all finite dimensions. A SIC- 
POVM is usually thought of as a highly symmetric structure in quantum state 
space. However, its elements can equally well be regarded as a basis for the Lie 
algebra g\(d, C). In this paper we examine the resulting structure constants, 
which are calculated from the traces of the triple products of the SIC-POVM 
elements and which, it turns out, characterize the SIC-POVM up to unitary 
equivalence. We show that the structure constants have numerous remarkable 
properties. In particular we show that the existence of a SIC-POVM in di- 
mension d is equivalent to the existence of a certain structure in the adjoint 
representation of gl(e?, C). We hope that transforming the problem in this way, 
from a question about quantum state space to a question about Lie algebras, 
may help to make the existence problem tractable. 
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1. Introduction 

Symmetric informationally complete positive operator- valued measures (SIC- 
POVMs) present us with what is, simultaneously, one of the most interesting, and 
one of the most difficult and tantalizing problems in quantum information [1-46]. 
SIC-POVMs are important practically, with applications to quantum tomography 
and cryptography [4, 8, 12, 15, 20, 29], and to classical signal processing [24, 36]. 
However, without in any way wishing to impugn the significance of the applications 
which have so far been proposed, it appears to us that the interest of SIC-POVMs 
stems less from these particular proposed uses than from rather broader, more gen- 
eral considerations: the sense one gets that SICs arc telling us something deep, 
and hitherto unsuspected about the structure of quantum state space. In spite of 
its being the central object about which the rest of quantum mechanics rotates, 
and notwithstanding the efforts of numerous investigators [47], the geometry of 
quantum state space continues to be surprisingly ill-understood. The hope which 
inspires our efforts is that a solution to the SIC problem will prove to be the key, 
not just to SIC-POVMs narrowly conceived, but to the geometry of state space in 
general. Such things are, by nature, unpredictable. However, it is not unreasonable 
to speculate that a better theoretical understanding of the geometry of quantum 
state space might have important practical consequences: not only the applica- 
tions listed above, but perhaps other applications which have yet to be conceived. 
On a more foundational level one may hope that it will lead to a much improved 
understanding of the conceptual message of quantum mechanics [7, 43, 45, 48]. 

Having said why we describe the problem as interesting, let us now explain why 
we describe it as tantalizing. The trouble is that, although there is an abundance of 
reasons for suspecting that SIC-POVMs exist in every finite dimension (exact and 
high-precision numerical examples [1, 2, 5, 11, 16, 19, 28, 39, 46] having now been 
constructed in every dimension up to 67) , and in spite of the intense efforts of many 
people [1-46] extending over a period of more than ten years, a general existence 
proof continues to elude us. In their seminal paper on the subject, published in 
2004, Renes et al [5] say "A rigorous proof of existence of SIC-POVMs in all finite 
dimensions seems tantalizingly close, yet remains somehow distant." They could 
have said the same if they were writing today. 

The purpose of this paper is to try to take our understanding of SIC mathematics 
(as it might be called) a little further forward. The research we report began with 
a chance numerical discovery made while we were working on a different problem. 
Pursuing that initial numerical hint we uncovered a rich and interesting set of 
connections between SIC-POVMs in dimension d and the Lie Algebra gl(d, C). The 
existence of these connections came as a surprise to us. However, in retrospect it 
is, perhaps, not so surprising. Interest in SIC-POVMs has, to date, focused on the 
fact that an arbitrary density matrix can be expanded in terms of a SIC-POVM. 
However, a SIC-POVM in dimension d does in fact provide a basis, not just for the 
space of density matrices, but for the space of all d x d complex matrices — i.e. the 
Lie algebra gl(d, C). Boykin et al [49] have recently shown that there is a connection 
between the existence problem for maximal sets of MUBs (mutually unbiased bases) 
and the theory of Lie algebras. Since SIC-POVMs share with MUBs the property 
of being highly symmetrical structures in quantum state space it might have been 
anticipated that there are also some interesting connections between SIC-POVMs 
and Lie algebras. 
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Our main result (proved in Sections 3, 4 and 5) is that the proposition, that a 
SIC-POVM exists in dimension d, is equivalent to a proposition about the adjoint 
representation of gl(d, C). Our hope is that transforming the problem in this way, 
from a question about quantum state space to a question about Lie algebras, may 
help to make the SIC-existence problem tractable. But even if this hope fails to 
materialize we feel that this result, along with the many other results we obtain, 
provides some additional insight into these structures. 

In d dimensional Hilbert space Oig a SIC-POVM is a set of d 2 operators E\, 
. . . of the form 

E r = in r (1) 
where the n r are rank-1 projectors with the property 

Tr(n r n s ) = | 1 1 r = S (2) 



d+1 



r 



We will refer to the U r as SIC projectors, and we will say that {II r : r — 1, . . . , d 2 } 
is a SIC set. 

It follows from this definition that the E r satisfy 



Y l E r =I (3) 



(so they constitute a POVM), and that they are linearly independent (so the POVM 
is informationally complete). 

It is an open question whether SIC-POVMs exist for all values of d. However, 
examples have been constructed analytically in dimensions 2-15 inclusive [1, 2, 11, 
16, 19, 28, 39, 46], and in dimensions 19, 24, 35 and 48 [16, 46]. Moreover, high 
precision numerical solutions have been constructed in dimensions 2-67 inclusive [5, 
46]. This lends some plausibility to the speculation that they exist in all dimensions. 
For a comprehensive account of the current state of knowledge in this regard, and 
many new results, see the recent study by Scott and Grassl [46]. 

All known SIC-POVMs have a group covariance property. In other words, there 
exists 

(1) a group G having d 2 elements 

(2) a projective unitary representation of G on Jf^: i.e. a map g — > U g from G 
to the set of unitaries such that U gi U g2 ~ U gig ,, for all <7i, <?2 (where the 
notation "~" means "equals up to a phase" ) 

(3) a normalized vector \ip) (the fiducial vector) 

such that the SIC-projectors are given by 

n g = u g W(ml ( 4 ) 

(where we label the projector by the group element g, rather than the integer r as 
above). 

Most known SIC-POVMs arc covariant under the action of the Wcyl-Heisenbcrg 
group (though not all — see Renes et al [5] and, for an explicit example of a non 
Weyl-Heisenberg SIC-POVM, Grassl [19]). Here the group is Z d x Z d , and the 
projective representation is p — > D p , where p = (^1,^2) £ x Z d and D p is the 
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corresponding Weyl-Heisenberg displacement operator 

d-1 

D p = J2r (2r+Pl)P2 \ r +Pi)( r \ ( 5 ) 

r 

In this expression r = e ~a the vectors |0), . . . \d— 1) are an orthonormal basis, 
and the addition in |r+pi) is modulo d. For more details see, for example, ref. [16]. 

One should not attach too much weight to the fact that all known SIC-POVMs 
have a group covariance property as this may only reflect the fact that group co- 
variant SIC-POVMs are much easier to construct. So in this paper we will try to 
prove as much as we can without assuming such a property. One potential benefit 
of this attitude is that, by accumulating enough facts about SIC-POVMs in general, 
we may eventually get to the point where we can answer the question, whether all 
SIC-POVMs actually do have a group covariance property. 

The fact that the d 2 operators n r are linearly independent means that they form 
a basis for the complex Lie algebra gl(d, C) (the set of all operators acting on "Kj)- 
Since the II r are Hermitian, then iIT r forms a basis also for the real Lie algebra 
u(d) (the set of all anti-Hermitian operators acting on JQ). So for any operator 
A G gl(d, C) there is a unique set of expansion coefficients a r such that 

d 2 

A = a, - n '- ( 6 ) 

r=l 

To find the expansion coefficients we can use the fact that 



d 2 



d+ 1 „ 1 



^ Tr(n r n s ) ^—6 st --)=S rt (7) 

S — 1 

from which it follows 



1 Tr(IL^) - - Tr(A) (8) 



d d 
Specializing to the case A = H r IL s we find 

where 

T rst = Tr(ILTI s IL) (10) 

To a large extent this paper consists in an exploration of the properties of these 
important quantities, which wc will refer to as the triple products. They are inti- 
mately related to the geometric phase, in which context they are usually referred 
to as 3- vertex Bargmann invariants (see Mukunda et al [50], and references cited 
therein). We have, as an immediate consequence of the definition, 

rji rr\ rr) rrv* rri* rri* / -f -i \ 

rst trs str rts tsr srt \ J 

It is convenient to define 

Jrst = ~i (T rst — T rst ) (12) 

Rrst = 1 (^rst + T rst ) (13) 
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So J rst is imaginary and completely anti-symmetric; R rs t is real and completely 
symmetric. Both these quantities play a significant role in the theory. It follows 
from Eq. (9) that 

d 2 

[n r ,n s ] = J2 Jr stilt (14) 
*=i 

So the J rs t are structure constants for the Lie algebra gl(d, C). As an immediate 
consequence of this they satisfy the Jacobi identity: 

d 2 

^^(JrsbJtba + JstbJrba + JtrbJsba) = (15) 
6=1 

for all r, s, £, a. The Jacobi identity holds for any representation of the structure 
constants. In the following sections we will derive many other identities which are 
specific to this particular representation. 

Turning to the quantities R rs t, it follows from Eq. (9) that they feature in the 
expression for the anti-commutator 

{n r , nj = J2 R rst u t - (is) 

They also play an important role in the description of quantum state space. Let 
p be any density matrix and let p r = \ Tr(H r p) be the probability of obtaining 
outcome r in the measurement described by the POVM with elements glL-. Then 
it follows from Eq. (8) that p can be reconstructed from the probabilities by 

d " ( 1\ 

P = ]Tf(d + iK--jii r (17) 

Suppose, now, that the p r are any set of d 2 real numbers. So we do not assume 
that the p r are even probabilities, let alone the probabilities coming from a density 
matrix according to the prescription p r = \ Tr(II r( o). Then it is shown in ref. [34] 
that the p r are in fact the probabilities coming from a pure state if and only if they 
satisfy the two conditions 



IX = 5(5+1) (18) 



r=l 

d 2 



£ RrstPrPsPt = (19) 
r,s,t—l v 1 

Let us look at the quantities J rs t and R rs t in a little more detail. For each r 
choose a unit vector \ip r ) such that II r = \ip r )(ip r \. Then the Gram matrix for these 
vectors is of the form 

G rs = (Vr# S ) = Krse"" (20) 

where the matrix 9 rs is anti-symmetric and 

' dS rs + 1 



< 21 > 

Note that the SIC-POVM does not determine the angles 9 rs uniquely since making 
the replacements \ip r ) — > e'^IVv) leaves the SIC-POVM unaltered, but changes 
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the angles 9 rs according to the prescription 9 rs — > 8 rs — (f> r + <p s . This freedom 
to rephase the vectors \ip r ) is not usually important. However, it sometimes has 
interesting consequences (see Section 9). It can be thought of as a kind of gauge 
freedom. 

The Gram matrix satisfies an important identity. Every SIC-POVM has the 
2-dcsign property [5, 17] 

d2 2d 

n r ® n r = J^Psym (22) 

r— 1 

where P sym is the projector onto the symmetric subspace of Jf^ <g> JCi. Expressed 
in terms of the Gram matrix this becomes 
d 2 d 

'^2G Sir G S2r G r t 1 G r t 2 = ^ — - ( y G Sl t 1 G S2 t 2 + G Sl t 2 Gsvti) (23) 

r— 1 

Turning to the triple products we have 

T rs t = G rs G stGtr — K rs K st K tr e ierst (24) 

where 

drst = @rs + ®st + ®tr (25) 

Note that the tensor 9 rst is completely anti-symmetric. In particular 9 rst — if any 
two of the indices are the same. Note also that re-phasing the vectors \ip r ) leaves 
the tensors T rst and 9 rst unchanged. They are in that sense gauge invariant. 
Finally, we have the following expressions for J rst and R rs t'- 

2i 

Jrst = -T^=T Sin Orat (26) 

d\>a + 1 

Rrst = —-T—K rs K at Ktr COB 9 ra t (27) 

Like the triple products, J rst and R rs t are gauge invariant. 

For later reference let us note that the matrix J r , with matrix elements 

{Jr)st = Jrst (28) 

is the adjoint representative of IT r in the SIC-projector basis: 

d 2 

adn r n s = [n r , n s ] = ^ j rst n t (29) 
t=i 

It can be seen that all the interesting features of the tensor G rs (respectively, 
the tensors T rstl J rs t and R rs t) are contained in the order-2 angle tensor 9 rs (re- 
spectively, the order-3 angle tensor 9 rst ). It is also easy to see that, for any unitary 
U, the transformation 

n. r -> J7n r c/ t (30) 

leaves the angle tensors invariant. This suggests that we shift our focus from indi- 
vidual SIC-POVMs to families of unitarily equivalent SIC-POVMs— SIC-families, 
as we will call them for short. 

Wc begin our investigation in Section 2 by giving necessary and sufficient con- 
ditions for an arbitrary tensor 9 rs (respectively 9 rs t) to be the rank-2 (respectively 
rank-3) angle tensor corresponding to a SIC-family. We also show that either angle 
tensor uniquely determines the corresponding SIC-family. Finally we describe a 
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method for reconstructing the SIC-family, starting from a knowledge of either of 
the two angle tensors. 

In Sections 3, 4 and 5 we prove the central result of this paper: namely, that 
the existence of a SIC-POVM in dimension d is equivalent to the existence of a 
certain very special set of matrices in the adjoint representation of gl(d, C). In 
Section 3 we show that, for any SIC-POVM, the adjoint matrices J r have the 
spectral decomposition 

J r = Qr - Qj (31) 

where Q r is a rank d — 1 projector which has the remarkable property of being 
orthogonal to its own transpose: 

Q r Qj = (32) 

We refer to this feature of the adjoint matrices as the Q-Q T property. In Section 3 
we also show that from a knowledge of the J matrices it is possible to reconstruct 
the corresponding SIC-family. In Section 4 we characterize the general class of 
projectors which have the property of being orthogonal to their own transpose. 
Then, in Section 5, we prove a converse of the result established in Section 3. The 
Q-Q T property is not completely equivalent to the property of being a SIC set. 
However, it turns out that it is, in a certain sense, very nearly equivalent. To be 
more specific: let L r be any set of d 2 Hermitian operators which constitute a basis 
for g\(d, C) and let C r be the adjoint representative of L r in this basis. Then the 
necessary and sufficient condition for the C r to have the spectral decomposition 

Cr = Qr - Ol (33) 

where Q r is a rank d — 1 projector such that Q r Qj = is that there exists a 
SIC set II,. such that L r = e r (n r + al) for some fixed number a € M and signs 
e r = ±1. In particular, the existence of an Hermitian basis for gl(d, C) having the 
Q-Q T property is both necesary and sufficient for the existence of a SIC-POVM in 
dimension d. 

In Section 6 we digress briefly, and consider s\(d, C) (the Lie algebra consisting 
of all trace-zero d x d complex matrices). As we have explained, this paper is 
motivated by the hope that a Lie algebraic perspective will cast light on the SIC- 
existence problem, rather than by an interest in Lie algebras as such. We focus on 
gl(d, C) because that is the case where the connection with SIC-POVMs seems most 
straightforward. However a SIC-POVM also gives rise to an interesting geometrical 
structure in sl(d, C), as we show in Section 6. 

In Section 7 we derive a number of additional identities satisfied by the J and 
Q matrices. 

The complex projectors Q r , Qj and the real projector Q r + Qj define three 
families of subspaces. It turns out that there are some interesting geometrical 
relationships between these subspaces, which we study in Section 8. 

Finally, in Section 9 we show that, with the appropriate choice of gauge, the 
Gram matrix corresponding to a Weyl-Heisenberg covariant SIC-family has a fea- 
ture analogous to the Q-Q T property, which we call the P-P T property. It is an 
open question whether this result generalizes to other SIC-families, not covariant 
with respect to the Weyl-Heisenberg group. 
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2. The Angle Tensors 

The purpose of this section is to establish the necessary and sufficient conditions 
for an arbitrary tensor 6 rs (respectively 9 rst ) to be the order-2 (respectively order- 
3) angle tensor for a SIC-family. We will also show that either one of the angle 
tensors is enough to uniquely determine the SIC-family. Moreover, we will describe 
explicit procedures for reconstructing the family, starting from a knowledge of one 
of the angle tensors. 

We begin by considering the general class of POVMs (not just SIC-POVMs) 
which consist of d 2 rank-1 elements. A POVM of this type is thus defined by a set 
of d 2 vectors ■ ■ • , |6 2 ) with the property 

d 2 

!>><£r|=i (34) 
r=l 

d 2 II 1 1 2 

Note that J^r=i 16} = ^> so ^ ne vec tors |£ r ) cannot all be normalized. In the 
particular case of a SIC-POVM the vectors all have the same norm |||6)|| = tt|- 
However in the general case they may have different norms. 
Given a set of such vectors consider the Gram matrix 

Pre = (616) (35) 

Clearly the Gram matrix cannot determine the POVM uniquely since if U is any 
unitary operator then the vectors U\£ r ) will define another POVM having the same 
Gram matrix. However, the theorem we now prove shows that this is the only free- 
dom. In other words, the Gram matrix fixes the POVM up to unitary equivalence. 
The theorem also provides us with a criterion for deciding whether an arbitrary 
d 2 x d 2 matrix P is the Gram matrix corresponding to a POVM of the specified 
type. As a corollary this will give us a criterion for deciding whether an arbitrary 
tensor 9 rs is specifically the ordcr-2 angle tensor for a SIC-family. 

Theorem 1. Let P be any d 2 x d 2 Hermitian matrix. Then the following conditions 
are equivalent: 

(1) P is a rank d projector. 

(2) P satisfies the trace identities 

Tr(P) = Tr(P 2 ) = Tr(P 3 ) = Tr(P 4 ) = d (36) 

(3) P is the Gram matrix for a set of d vectors |6) (not all normalized) such 
that |6) (61 is a POVM: 

(616) =Prs (37) 



£16X61 = / 



(38) 



Suppose P satisfies these conditions. To construct a POVM corresponding to P 
let the d column vectors 



(39) 













62 


1 


62 


1 • • • 1 


62 










\€dd 2 J 
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be any orthonormal basis for the subspace onto which P projects. Define 

d 

IU=£C|o) (40) 

a=l 

where the vectors \a) are any orthonormal basis for JLj. Then P is the Gram matrix 
for the vectors ■ ■ ■ > \£d 2 )- Moreover, the necessary and sufficient condition for 
any other set of vectors \r]i), . . . , \rj,p) to have Gram matrix P is that there exist a 
unitary operator U such that 

\Vr) = U\£ r ) (41) 

for all r. 

Proof. We begin by showing that (3) => (1). Suppose • • ■ \£,d 2 ) is an Y set of 
d 2 vectors such that |£ r )(£r| is a POVM. So 



El^r>«r| = I 



r=l 



Let 



be the Gram matrix. Then P is Hermitian. Moreover, P 2 = P since 



(42) 
(43) 



1 = 1 



P 



Also 



Tr(p) = E(erie r > = d 



(44) 



(45) 



(as can be seen by taking the trace on both sides of Eq. (42)). So P is a rank-c? 
projector. 

We next show that (1) => (3). Let P be a rank-d projector, and let the d 
column vectors 













62 


1 


£.22 


, . . . , 


S,d2 






{.fa) 




\£dd 2 / 



(46) 

be an orthonormal basis for the subspace onto which it projects. So 

d 2 

EC4r = <U (47) 

r=l 

for all a, 6, and 

d 2 

Y,^rCs=Prs (48) 
a=l 



9 



for all r, s. Now let . . . |£<p) be the vectors defined by Eq. (40). Then it follows 
from Eq. (47) that 



r=l 



Eifr>«ri= E Ece 

a, 6=1 \l 

i>>< 



*)<6| 



a,6=l \r=l 
d 

J 



implying that |£ r )(£ r | is POVM. Also, it follows from Eq. (48) that 



\SrlSs/ / ^ Sar>as r. 



(49) 



(50) 



implying that the |£ r ) have Gram matrix P. 

We next turn to condition (2). The fact that (1) =$> (2) is immediate. To 
prove the reverse implication observe that condition (2) implies 

Tr(P 4 ) - 2 Tr(P 3 ) + Tr(P 2 ) = (51) 

Let Ai, . . . , Ad 2 be the eigenvalues of P. Then Eq. (51) implies 







(52) 



It follows that each eigenvalue is either or 1. Since Tr(P) = d we must have d 
eigenvalues = 1 and the rest all = 0. So P is a rank-d projector. 

It remains to show that the POVM corresponding to a given rank-d projector 
is unique up to unitary equivalence. To prove this let P be a rank-d projector, let 
|£ r ) be the vectors defined by Eq. (40), and let |r/i), . . . , \r]d 2 ) be any other set of 
vectors such that 

(Vr\Vs) = P rs 



for all r, s. Define 
Then 



Var = (Vr\a) 



E^%r = ( a \ E MiVrl \b) = §ab 



(53) 
(54) 

(55) 



(because \rj r )(r] r \ is a POVM) and 



^ ] VarVas — ^rs 
a=l 

(because the \rj r ) have Gram matrix P). So the d column vectors 



fmi\ 

r)i2 



( 7/21 \ 



Vd2 



(56) 



(57) 
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are an orthonormal basis for the subspace onto which P projects. But the column 
vectors 



(58) 



are also an orthonormal basis for this subspace. So there must exist a dx d unitary 
matrix U a b such that 

d 







£.22 


, . . . , 


Cd2 






{.fa) 




\idd 2 J 



T]ar ^ ^ ^ab£br 



6=1 



for all a,r. Define 



Then 



U= U : b \a)(b\ 

a, 6=1 
\Vr) = U\£ r ) 



for all r. 

In the case of a SIC-POVM we have 



(59) 



(60) 



(61) 
□ 



(62) 



where the vectors \ip r ) are normalized, and 

Prs = \G rs = -K rs e ie " (63) 
d d 

where G is the Gram matrix of the vectors \tp r ) and 9 rs is the order-2 angle tensor. 
In the sequel we will distinguish these matrices by referring to G as the Gram 
matrix and P as the Gram projector. 
We have 

Corollary 2. Let 8 rs be a real anti- symmetric tensor. Then the following state- 
ments are equivalent: 

(1) 8 rs is an order-2 angle tensor corresponding to a SIC-family. 

(2) 9 rs satisfies 

d 2 

J2 K rt K ts e^+ 8 ^ = dK rs e w " (64) 



for all r, s. 
(3) 9 rs satisfies 



J2 K rs K st K tr e^+ e ^) = d * 



(65) 



r.s.t— 1 



r,s,t,u—l 



(66) 
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Let H r , Tl' r be two different SIC-sets, and let O rs , 9' rs be corresponding order-2 
angle tensors. Then there exists a unitary U such that 

n' r = UU r U^ (67) 

for all r if and only if 

e' rs = e rs - 4>r + <Ps (68) 

for some arbitrary set of phase angles <f> r (in other words two SIC-sets are unitarily 
equivalent if and only if their order-2 angle tensors are gauge equivalent). 

A SIC-family can be reconstructed from its order-2 angle tensor 9 rs by calculating 
an orthonormal basis for the subspace onto which the Gram projector 

Prs = \Kt,<P" (69) 

a 

projects, as described in Theorem 1. 

Remark. The sense in which we are using the term "gauge equivalence" is explained 
in the passage immediately following Eq. (21). 

Note that condition (2) imposes d 2 (d 2 — l)/2 independent constraints (taking 
account of the anti-symmetry of 6 rs ). Condition (3), by contrast, only imposes 2 
independent constraints. It is to be observed, however, that the price we pay for 
the reduction in the number of equations is that Eqs. (65) and (65) are respectively 
cubic and quartic in the phases, whereas Eq. (64) is only quadratic. 

Proof. Let 9 rs be an arbitrary anti-symmetric tensor, and define 

Prs = -J< rs e w ™ (70) 
a 

The anti-symmetry of 9 rs means that P is automatically Hermitian. So it follows 
from Theorem 1 that a necessary and sufficient condition for P rs to be a rank-d 
projector, and for 9 rs to be the order-2 angle tensor of a SIC-family, is that 

d 2 

K rt K ts e^+ e ^ = dK rs e w ™ (71) 

t=i 

for all r, s. 

To prove the equivalence of conditions (1) and (3) note that the conditions 
Tr(P) = Tr(P 2 ) = d are an automatic consequence of P having the specified form. 
So it follows from Theorem 1 that 9 rs is the order-2 angle tensor of a SIC-family if 
and only if Eqs. (65) and (66) are satisfied. 

Now let il r , W r be two SIC-sets and let 9 rs , 9' rs be ordcr-2 angle tensors corre- 
sponding to them. Then there exist normalized vectors \ip r ), \ip' r ) such that 

n r H^)(Vvl K = mm (72) 

for all r, and 

WM)=K r .e i0 '" (73) 
for all r, s. 

Suppose, first of all, that there exists a unitary U such that 

n; = uu r u^ (74) 



(ip r \ip s ) = K rs e 
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Then there exist phase angles <p r such that 

\f// r ) = et'Utyr) (75) 
for all r, which is easily seen to imply that 



y rs ra tt i 



for all r, s. So 9 rs , 9' rs are gauge equivalent. 

Conversely, suppose there exist phase angles cj> r such that 



(76) 



(77) 



Dehnc 

K)=e-^|<) (78) 

Then 

WW:)=K rs e ie - = (i> r W s ) (79) 
for all r, s. So it follows from Theorem 1 that there exists a unitary U such that 

Wr) = U\i> r ) (80) 

for all r. Consequently 

n; = K)«| = [/n r [7t (8 i) 

for all r. So il r and Ilj, are unitarily equivalent. □ 

We now turn to the order-3 angle tensors. We have 

Theorem 3. Let 6 rst be a real completely anti- symmetric tensor. Then the follow- 
ing conditions are equivalent: 

(1) 9 rst is the order-3 angle tensor for a SIC-family 

(2) For some fixed a and all r, s,t 



7 ars ~r "ast t "atr — "rst 



(82) 



and for all r, s 

d 2 

K rt K ts e id <->* = dK rs (83) 

t=i 

(3) For some fixed a and all r,s,t 

Gars + &ast + datr = &rst (84) 

and 

d 2 

K rs K st K tr e w ^ = d 4 (85) 

r,s,t— 1 

d 2 

K rs K st K tu K ur e l ^ +Bt ^ = d 5 (86) 

r.s.t.u— 1 
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Let II r , H' r be two different SIC-sets and let 9 rst , 9' rst be the corresponding order- 
3 angle tensors. Then the necessary and sufficient condition for there to exist a 
unitary U such that 

n; = un r tf (87) 

for all r is that 9' rst = 9 rst for all r,s,t (in other words two SIC-sets are unitarily 
equivalent if and only if their order-3 angle tensors are identical). 

Let 9 rst be the order-3 angle tensor corresponding to a SIC-family. Then the 
order-2 angle tensor is given by (up to gauge freedom) 



for any fixed a, from which the SIC-family can be reconstructed using the method 
described in Theorem 1. 

Remark. Unlike the order-2 tensor, the order-3 angle tensor is gauge invariant. This 
means that it provides what is, in many ways, a more useful characterization of 
the SIC-family. For that reason we will be almost exclusively concerned with the 
ordcr-3 tensor in the remainder of this paper. 

Proof. The fact that (1) (2) is an immediate consequence of the definition of 

the order-3 angle tensor and condition (2) of Corollary 2. To prove that (2) (1) 
let 9 rs t be a completely anti-symmetric tensor such that condition (2) holds. Define 

9 rs = 9 ars (89) 

for all r, s. Then Eq. (83) implies 

d 2 / d 2 \ ' 

£ K rt K ts e^+^) = e «~ £ K rt K ts e^ = dK rs e if> ™ (90) 
t=i \ t=1 / 

for all r, s. It follows from Corollary 2 that 9 rs is the order-2 and 9 rst the order-3 
angle tensor of a SIC-family. 

The equivalence of conditions (1) and (3) is proved similarly. 

It remains to show that two SIC-sets are unitarily equivalent if and only if 
their order-3 angle tensors are identical. To see this let II r = \ip r )(ip r \ an d nj. = 
\ifj' r ) (ip' r \ be two different SIC-sets having the same order-3 angle tensor 9 rst . Let 
9 rs (respectively 9' rs ) be the order-2 angle tensor corresponding to the vectors \tp r ) 
(respectively \ip' r ))- Choose some fixed index a. We have 



for all r, s. Consequently 



for all r, s, where 



= Kr + Osa + Ore (91) 

= 9 rs + <t> r - <t> s (92) 
= e ar - 9' ar (93) 



So 9' rs and 9 rs are gauge equivalent. It follows from Corollary 2 that II,, and IIJ, are 
unitarily equivalent. Conversely, suppose that II,. and IL^, are unitarily equivalent, 
and let 9 rs , 9' rs be order-2 angle tensors corresponding to them. It follows from 
Corollary 2 that 9 rs and 9' rs are gauge equivalent. It is then immediate that the 
order-3 angle tensors are identical. □ 
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Finally, let us note that when expressed in terms of the triple products Eq. (83) 
reads 



J2T rst = dKl (94) 
t=i 

while Eq. (85) reads 

d 2 

]T T rst = d A (95) 

r,s 7 t — l 

For Eq. (86) we have to work a little harder. We have 



1 K r 

r,8 J t,u=± 



J2 -^r TrstTtur = d5 (96) 

r, 

from which it follows 

(ft = (-d$rt + d + l)T rst T tur 

r.s.t.u—l 



(d+l) J! T rst T tur -d J2 K rs K ru 
r,s,t,u—l r,s y u=l 

d 2 

(d+l) T rst T tur -d 5 (97) 



r,s,t 1 



Consequently 



d d 2d 5 

^ Tr(T s T„) = ^ T rst T tur = — — (98) 



d+l 

s,ii— 1 r,s,t.u— 1 



This equation be alternatively written 



d + i 

r,s—l 



where T r is the matrix with matrix elements (T r ) uv = T ruv . 

When they are written like this, in terms of the triple products, the fact that 
Eq. (94) implies Eqs. (95) and (98) becomes almost obvious. The reverse implica- 
tion, by contrast, is rather less obvious. 



3. Spectral Decompositions 
Let T r , J r , R r be the d 2 x d 2 matrices whose matrix elements are 

{T r ) s t = T rst (Jr)st = Jrst (Rr)st = Rrst (100) 

where J rs t, R rs t are the quantities defined by Eqs. (12) and (13). So J r is the adjoint 
representation matrix of n r . In this section we derive the spectral decompositions 
of these matrices. To avoid confusion we will use the notation \ip) to denote a ket in 
d dimensional Hilbert space Jf^, and to denote a ket in d 2 dimensional Hilbert 
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space "KcP ■ In terms of this notation the spectral decompositions will turn out to 
be: 

T r = -^Q r + ^L\\e r ))((e r \\ (101) 

J r = Q r - QJ (102) 

R r = Q r + Qj + A\\e r ))((e r \\ (103) 

In these expressions the vector \\e r )) is normalized, and its components in the stan- 
dard basis are all real. Q r is a rank d — 1 projector such that 

Q r \\e r ))=Q?\\e r ))=0 (104) 

and which has, in addition, the remarkable property of being orthogonal to its own 
transpose (also a rank d — 1 projector): 

Q r Qj = (105) 

Explicit expressions for ||e r )) and Q r will be given below. 
It will be convenient to define the rank 2(d — 1) projector 

Rr = Q r + QT (106) 

We have 

R r = J 2 (107) 

and 

J2 r = J R P + 4||e P ))<(e P || (108) 
Since Q r is Hcrmitian we have 

Ql = Q*r (109) 
where Q* is the matrix whose elements are the complex conjugates of the corre- 
sponding elements of Q r - So R r is twice the real part of Q r and — iJ r is twice the 
imaginary part. 

In Section 5 we will show that Eq. (102) is essentially definitive of a SIC-POVM. 
To be more specific, let L r be any set of d 2 Hcrmitian matrices which constitute a 
basis for gl(d, C), and let C r be the adjoint representative of L r in that basis. Then 
we will show that C r has the spectral decomposition 

a,, = Q r -Qj (no) 

where Q r is a rank d — 1 projector which is orthogonal to its own transpose if and 
only if the L r are a family of SIC projectors up to multiplication by a sign and 
shifting by a multiple of the identity. 

Having stated our results let us now turn to the task of proving them. We begin 
by deriving the spectral decomposition of T r . Multiplying both sides of the equation 

1 d ' 2 

-£r rst n t - K 2 rs i (in) 



d 
t=i 



by n r we find 

d 



n,,n s = ^±1 J2 r rst n r n t - K 2 rs n r 



{d+lf 
d 2 



We have 



1 d 

— -y^T rst (di5 rf + 1) 
1 / d2 

— I dT rsr + ^ T r. 



-T 



Consequently 



2d 
d+l' 

2d 
d + l~ 



Comparing with Eq. (Ill) we deduce 



{Trf rs - d + 1 T rst + d + 1 K^Jift 



Now define 



HO) 



d + l 
2d 



where the basis kets \\s)) are given by (in column vector form) 





f 1 ) 

















1 







lll» = 




,l|2» = 




,...,lld 2 )) = 






W 




w 




w 



It is easily verified that ||e r )) is normalized. Eq. (115) then becomes 

d „ 2d 2 



-T 



d+l r (d + l) 2 



2 ll^/A^rl 



Using Eq. (113) we find 



(s||T r ||e r )) = 



d+l 
2d 



t=l 
d+l TS 



2d 
d+l 



((s\\e r )) 



So j|e, )) is an eigenvector of T r with eigenvalue 
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Also define 



Q r = l±lr r -2|| ej .))((e r || (120) 
a 



So in terms of the order-3 angle tensor the matrix elements of Q r are 

d+1 



a 

Q r is Hermitian (because T r is Hermitian). Moreover 



K rs K rt (K st e lB -t - K rs K rt ) (121) 



Qr = ^^Tr - 4 He.))((er|| = Qr (122) 

So Q r is a projection operator. Since 

d 2 

Tr(T r ) - Yl T ™« = E = d (123) 

u u— 1 

we have 

Tr(Q r ) = d - 1 (124) 
We have thus proved that the spectral decomposition of T r is 

Tr = ^ ri Qr + ^- l \\er)){{e r \\ (125) 

where <5 r is a rank d — 1 projector, as claimed. 

We next prove that Qj\\e r )) = 0. The fact that the components of ||e r )) in the 
standard basis are all real means 

«,s||r r T || er )> - «e r ||T r || S » = -^{{s\\e r )) (126) 

So || e r )) is an eigenvector of Tj as well as T r , again with the eigenvalue ^pj. In 
view of Eq. (120) it follows that Qj\\e r )) = 0. 

Turning to the problem of showing that Q r is orthogonal to its own transpose. 
We have 



; r ^ r 



d 

d 2 

It follows from Eq. (24) that 



V-2IMKH) ( d -±±Tj -2\\e r )){{e r \ 



T r Tj -A\\e r )){{e r \\ (127) 



((s||T r r r T ||t)) = ^r rs j rfu 

u=l 

d 2 

= G rs G rt G su Gt u G ur G ur (128) 

u=l 

In view of Eq. (23) (i.e. the fact that every SIC-POVM is a 2-design) this implies 

({ S \\T r T?\\t)) = ^L\Grs\ 2 \G rt \ 2 

d+1 rs rt 
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So 



and consequently 



^^(( S ||e r ))((e r ||t)) (129) 



T r T? = j J —^\\e r ))((e r \\ (130) 



Q r Qj = (131) 



Eqs. (102) and (103) are immediate consequences of the results already proved 
and the definitions of J r , R r . 

We defined the J matrices to be the adjoint representatives of the SIC-projectors, 
considered as a basis for the Lie algebra gl(d, C), and that is certainly a most 
important fact about them. However, the results of this section show that, along 
with the vectors ||e r )), they actually determine the whole structure. Specifically, 
we have 

Qr=\{Jr + Jr) (132) 

Rr = J ? 2 +4|| er ))((e r || (133) 

T r = ^pyy (J r + J 2 r + 4||e r )> «e r H) (134) 

Moreover, if we know the T matrices then we know the order-3 angle tensor, which 
in view of Theorem 3 means we can reconstruct the SIC-projectors. Since the 
vectors ||e r )) are given, once and for all, this means that the problem of proving the 
existence of a SIC-POVM in dimension d is equivalent to the problem of proving the 
existence of a certain remarkable structure in the adjoint representation of gl(d, C) 
(as we will see in more detail in Section 5). 

In the Introduction we began with the concept of a SIC-POVM, and then defined 
the J matrices in terms of it. However, one could, if one wished, go in the opposite 
direction, and take the Lie algebraic structure to be primary, with the SIC-POVM 
being the secondary, derivative entity. 



4. The Q-Q t Property 

The next five sections are devoted to a study of the J matrices which, as we will 
see, have numerous interesting properties. We begin our investigation by trying to 
get some additional insight into what we will call the Q-Q T property: namely, the 
fact that the J matrices have the spectral decomposition 

J r = Q r -Qj (135) 

where Q r is a rank d— 1 projector which is orthogonal to its own transpose. We wish 
to characterize the general class of matrices which arc of this type. The following 
theorem provides one such characterization. 

Theorem 4. Let A be a Hermitian matrix. Then the following statements are 
equivalent: 
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(1) A has the spectral decomposition 

A = P-P T 

where P is a projector which is orthogonal to its own transpose. 

(2) A is pure imaginary and A 2 is a projector. 



(136) 



Proof. To show that (1) =>■ (2) observe that the fact that P is Hcrmitian means 

where P* is the matrix whose elements are the complex conjugates of the corre- 
sponding elements of P. So Eq. (136) implies that the components of A are pure 
imaginary. Since PP T = it also implies that A 2 is a projector. 

To show that (2) =>■ (1) observe that the fact that A 2 is a projector means 
that the eigenvalues of A = ±1 or 0. So 

A = P-P' (138) 

where P, P' arc orthogonal projectors. Since A is pure imaginary we must have 

p t _ (p/yr = A T = A* = -A = P' - P (139) 

P T and (P') T are also orthogonal projectors. So if P T \ip) = \ip), and \4>) is nor- 
malized, we must have 

1 = (V|P T |V) 

= <V>|(P T -(P') T )|V) 

= (VIP'I^) - (V>|P|V>) (140) 

Since 

< (V>|P'|V>> < 1 (141) 
< {il)\P\ip) < 1 (142) 

we must have (ifi\P'\ip) = 1, implying = IV')- Similarly P'\ip) = implies 

P T |V) = |V>. So 

P' = P T (143) 

□ 

We also have the following statement, inspired in part by Ref. [51], 

Theorem 5. The necessary and sufficient condition for a matrix P to be a projector 
which is orthogonal to its own transpose is that 

P = SDS T (144) 

where S is an any real orthogonal matrix and D has the block- diagonal form 

(a ... ... 0\ 



D = 








a 











(145) 
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with 



1 [1 -i 
1 



(146) 



In other words D has n copies of a on the diagonal, where n = rank(P), and 
everywhere else. 

Proof. Sufficiency is an immediate consequence of the fact that a is a rank 1 pro- 
jector such that acr T = 0. 

To prove necessity let d be the dimension of the space and n the rank of P. It 
will be convenient to define 



|1) = 







W 



12} = 



1 



\d) 







Vv 



(147) 



In terms of these basis vectors we have 



P= J2 P rs\r)(s\ 



(148) 



Now let |ai), . . . , |o n ) be an orthonormal basis for the subspace onto which P 
projects, and let \a*) be the column vector which is obtained from \a r ) by tak- 
ing the complex conjugate of each of its components. Taking complex conjugates 
on each side of the equation 



gives 



So 



P\a r ) = \a r ) 



P*\a* 



(149) 



(150) 



, | a*) is an orthonormal basis for the subspace onto which P T = P* 



projects. Since P is orthogonal to P we conclude that 



for all r, s. 

Next define vectors 



(a r K)=0 
• ■ , \b2n) by 

|6 2r -i> = ^§(K>-k-» 

l 6 2r> = 4?(K> + K» 



(151) 

(152) 
(153) 



By construction these vectors are orthonormal and real. So we can extend them 
to an orthonormal basis for the full space by adding a further d — In vectors 
\b2n+i), ■ ■ ■ j \bd)> which can also be chosen to be real. We have 

n 

P = |a r )(a r | 

r=l 
1 " 

= x$^(|&ar-i)<&ar-i| -i\b2r-x)(h r \ +i\b 2r ){b 2r - x \ + \b 2r ){b 2r \) (154) 



r=l 
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So if wc define 



S = £»<7 



(155) 



then S is a real orthogonal matrix such that 



P = SDS 1 



(156) 



where 



1 " 

D = -^(l 2r ~ 1 ^ 2r ~ 1 ' ~*' 2r ~ 1 ^ 2r ' + i l 2r )( 2r - !! + l 2r, )( 2r l) ( 157 ) 



is the matrix defined by Eq. (145). 



□ 



This result implies the following alternative characterization of the class of ma- 
trices to which the J matrices belong 

Corollary 6. Let A be a Hermitian matrix. Then the following statements are 
equivalent: 

(1) A has the spectral decomposition 

A = P-P T (158) 

where P is a projector which is orthogonal to its own transpose. 

(2) There exists a real orthogonal matrix S such that 



A = SDS T 
where D has the block diagonal form 



(159) 



D = 






Vo 

<jy being the Pauli matrix 







<Jy 









(160) 



0/ 







(161) 



In other words D has n copies of o~ y on the diagonal, where n = | rank(A) , 
andO everywhere else (note that a matrix of this type must have even rank). 



Proof. Immediate consequence of Theorem 5. 

5. Lie Algebraic Formulation of the Existence Problem 



□ 



This section is the core of the paper. We show that the problem of proving the 
existence of a SIC-POVM in dimension d is equivalent to the problem of proving 
the existence of an Hermitian basis for g\(d, C) all of whose elements have the Q-Q T 
property. Wc hope that this new way of thinking will help make the SIC-existence 
problem more amenable to solution. 

The result we prove is the following: 
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Theorem 7. Let L r be a set of d 2 Hermitian matrices forming a basis for gl(d, C). 
Let C rst be the structure constants relative to this basis, so that 

d 2 

[L r ,L s }=J2 C ™tLt (162) 

i=l 

and let C r be the matrix with matrix elements (C r ) st = C rst . Then the following 
statements are equivalent 

(1) Each C r has the spectral decomposition 

C r = P r - Pj (163) 

where P r is a rank d—1 projector which is orthogonal to its own transpose. 

(2) There exists a SIC-set IT r , a set of signs e r = ±1 and a real constant 
a 7^ — \ such that 

L r = e r (n r + aL) (164) 

Remark. The restriction to values of a ^ — \ is needed to ensure that the matrices 
L r are linearly independent, and therefore constitute a basis for gl(d, C) (otherwise 
they would all have trace = 0). The Q-Q T property continues to hold even if a 
does = — 2- 

It will be seen that it is not only SIC-sets which have the Q-Q T property, but 
also any set of operators obtained from a SIC-set by shifting by a constant and 
multiplying by an r-dcpendent sign. So the Q-Q T property is not strictly equivalent 
to the property of being a SIC-sct. However, it could be said that the properties 
are almost equivalent. In particular, the existence of an Hermitian basis for gl(e£, C) 
having the Q-Q T property implies the existence of a SIC-POVM in dimension d, 
and conversely. 

Proof that (2) => (1). Taking the trace on both sides of 

d 2 

[n r ,n s ] =jT j rst n t (165) 
t=i 

we deduce that 

d 2 

J rst = (166) 

t=l 

Then from the definition of L r in terms of II r we find 

Crst = ^r e s^tJrst (167) 

Consequently 

C r = P r - Pj (168) 

where 

P r = e r SQ r S (169) 
S being the symmetric orthogonal matrix 

/ei . . . \ 
e 2 ... 



S 



(170) 



\0 ... £ d 2) 



The claim is now immediate. 
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Proof that (1) =>■ (2). For this we need to work harder. Since the proof is rather 
lengthy we will break it into a number of lemmas. We first collect a few elementary 
facts which will be needed in the sequel: 

Lemma 8. Let L r be any Hermitian basis for g\(d, C), and let C rst and C r be 
the structure constants and adjoint representatives as defined in the statement of 
Theorem 7. Let l r = Tr(L r ). Then 

(1) The l r are not all zero. 

(2) The C rs t are pure imaginary and antisymmetric in the first pair of indices. 

(3) The C rst are completely antisymmetric if and only if the C r are Hermitian. 

(4) Ln every case 



]T C rst l t = (171) 




t=l 

for all r, s. 

(5) In the special case that the C r are Hermitian 
d 2 

J2l r L r = Kl (172) 
r=l 

where 

1 \ 

(173) 

Proof. To prove (1) observe that if the l r were all zero it would mean that the 
identity was not in the span of the L r — contrary to the assumption that they form 
a basis. 

To prove (2) observe that taking Hermitian conjugates on both sides of Eq. (162) 
gives 

d 2 

— [L r , L s ] = ^2 C* st L t (174) 
t=i 

from which it follows that C* st = —C rst . The fact that C srt = —C rst is an imme- 
diate consequence of the definition. 

(3) is now immediate. 

(4) is proved in the same way as Eq. (166). 

To prove (5) observe that if the C r arc Hermitian it follows from (2) and (3) that 



J2 l r G rst = (175) 



r=l 

for all s, t. Consequently the matrix 



]T l r L r (176) 
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commutes with everything. But the only matrices for which that is true arc multi- 
ples of the identity. It follows that 

d 2 

k-Lr = Kl (177) 

r=l 

for some real k. Taking the trace on both sides of this equation we deduce 

d 2 

^2l 2 r =dn (178) 

r=l 

The fact that k > is a consequence of this and statement (1). □ 

We next observe that if the C r have the Q-Q T property they must, in particular, 
be Hermitian. It turns out that that is, by itself, already a very strong constraint. 

Before stating the result it may be helpful if we explain the essential idea on 
which it depends. Although we have not done so before, and will not do so again, it 
will be convenient to make use of the covariant /contravariant index notation which 
is often used to describe the structure constants. Define the metric tensor 

M rs = Tr(L r L s ) (179) 

and let M rs be its inverse. So 

d 2 f 1 = 

Y,M rt M ts =M\ = { (180) 

We can use these tensors to raise and lower indices (we use the Hilbert-Schmidt 
inner product for this purpose because the fact that gl(d, C) is not semi-simple 
means that its Killing form is degenerate [52-55]). In particular, the matrices 

d 2 

L r = ]T M rs L s (181) 
t=i 

are the basis dual to the L r : 

Tr{L r L s ) = M r s (182) 

Suppose we now define structure constants C rs t by 

d 2 

[L r ,L s \=Y J C rs tL t (183) 
t=i 

(so in terms of the C rs t we have C rs t = C rst ). It follows from the relation 

C rst = Tr([L r ,L s ]L t ) = Tr(L r [L s ,L t ]) (184) 

that the C rs t are completely antisymmetric for any choice of the L r . If we now 
require that the matrices C r be Hermitian it means that, not only the C rst , but 
also the C rs t must be completely antisymmetric. Since the two quantities are 
related by 

d 2 

Crst = V"] C rS uM ut (185) 
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this is a very strong requirement. It means that the L r must, in a certain sense, 
be close to orthonormal (relative to the Hilbert-Schmidt inner product). More 
precisely, it means we have the following lemma: 

Lemma 9. Let L r , C rs t and C r be defined as in the statement of Theorem 7, and 
let l r = Tr(L r ). Then the C r are Hermitian if and only if 

Tr(L r L s ) = I35 rs + -yl r l s (186) 

where j3, 7 are real constants such that f3 > and 7 < \ ■ 
If this condition is satisfied we also have 



r=l 

d 



E«-I^F <188) 



r=l 

Proof. To prove sufficiency observe that, in view of Eq. (185), the condition means 

d 2 

Crst = PC rst + ~fk ^ C rsu l u (189) 
u=l 

In view of Lemma 8, and the fact that (3 7^ 0, this implies 

Crst = -nCrst (190) 

Since the C rs t are completely antisymmetric we conclude that the C rs t must be 
also. It follows that the C r are Hermitian. 

To prove necessity let C r (respectively M) be the matrix whose matrix elements 
are C rs t (respectively M st ). Then Eq. (185) can be written 

C r = C r M (191) 

Taking the transpose (or, equivalently, the Hermitian conjugate) on both sides of 
this equation we find 

C r = MC r (192) 

implying 

[M, C r ] = (193) 
for all r. Since the L r are a basis for g\(d, C) we deduce 

[M,ad A ]=0 (194) 

for all A € gl(d, C). Eq. (186) is a straightforward consequence of this, the fact 
that gl(d, C) has the direct sum decomposition CI © sl(d, C), the fact that sl(d, C) 
is simple, and Schur's lemma [52-55]. However, for the benefit of the reader who is 
not so familiar with the theory of Lie algebras we will give the argument in a little 
more detail. 
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Given arbitrary A = X)r=i a rL r , let \\A)) denote the column vector 



So 



\L r 







In view of Lemma 8 we then have 



1 



\L d 2) 



M 





(195) 



(196) 



(197) 



Since 



Tr(A) = ^a r l r = 



we have that A £ sl(d, C) if and only if ((I\\A)) = 0. 

Now observe that it follows from Lemma 8 and the definition of M that 

M\\I)) = n\\I)) 



(198) 



(199) 



If M is a multiple of the identity we have M rs = K,6 rs and the lemma is proved. 
Otherwise M has at least one more eigenvalue, j3 say. Let £ be the corresponding 
eigenspace. Since £ is orthogonal to ||/)) it follows from Eq. (198) that £ C sl(d, C). 
Since M commutes with every adjoint representation matrix we have 



adyi£ C £ 



(200) 



for all A £ s\(d, C). So £ is an ideal of sl(d, C). However sl(d, C) is a simple Lie 
algebra, meaning it has no proper ideals [52-55]. So we must have £ = s\(d, C). It 
follows that if we define 



then 



d 



M\\L r )) = ±M\\I)) + M\\L r )) 



d 
kI 



(201) 

(202) 
(203) 



Y^iP&ra + llrQ \\L S 



(204) 



s=l 



where 7 = \ (l - f J. Eqs. (186), (187) and (188) are now immediate (in view of 
Lemma 8). 
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It remains to establish the bounds on j3, 7. Let A = ^r=i a rL r be any non-zero 

d 2 

element of sl(d, C). Then J2 r =i a ^ r = ^' so m v i ew °f Eq. (186) we have 

d 2 

< Tr(A 2 ) = ft'^2 a l (205) 

It follows that j3 > 0. Also, using Lemma 8 once more, we find 

f * 

ir = ~y^Js Ti(L r L s ) 

K s=l 

K K ^ S 

s=l 

= l r f^ + djj (206) 
Since the l r cannot all be zero this implies 



K 



= 1 - *y (207) 



Since ^ > we deduce that 7 < 3- □ 

Eq. (186) only depends on the C r being Hermitian. If we make the assumption 
that the C r have the Q-Q T property we get a stronger statement: 

Corollary 10. Let L r , C rs t and C r be as defined in the statement of Theorem 7. 
Suppose that the C r have the spectral decomposition 

C r = P r - Pj (208) 

where P r is a rank d — 1 projector which is orthogonal to its own transpose. Then 

(1) For all r 

Tr{L r ) = e' r l (209) 

(2) For all r, s 

Tr(L^) = ^, + ^--^) (210) 

(3) 

d 2 

^e' r L r = dlI (211) 

r=l 

for some real constant I > and signs e' r = ±1. 

Proof. The proof relies on the fact that the Killing form for gl(rf, C) is related to 
the Hilbert-Schmidt inner product by [55] 

Tr(ad A adB) = 2dTr(AB) - 2Tr(A) Tr(B) (212) 

Specializing to the case A = B = L r and making use of the Q-Q T property we find 

d- 1 = dTr(L%) -ll (213) 

Using Lemma 9 we deduce 

ll = d l^±l (214 ) 
1 — d-f 
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It follows that 

l r = e' r l (215) 
for some real constant I > and signs e' r = ±1. The fact that the L r are a basis 
for g\(d, C) means the l r cannot all be zero. So we must in fact have I > 0. Using 
this result in Eq. (188) we find 

+ d 2 l 2 j = dl 2 (216) 

while Eq. (214) implies 

d(3 + dl 2 1 = d-l + l 2 (217) 

This gives us a pair of simultaneous equations in j3 and 7. Solving them we obtain 

'"ITI < 218 > 

-=^( i2 -jii) < 219 > 

Substituting these expressions into Eqs. (186) and (187) we deduce Eqs. (210) 
and (211). □ 

The next lemma shows that each L r is a linear combination of a rank-1 projector 
and the identity: 

Lemma 11. Let L be any Hermitian matrix £ gl(d, C) which is not a multiple of 
the identity. Then 

rank(ad L ) > 2(d - 1) (220) 

The lower bound is achieved if and only if L is of the form 

L = nI + ^P (221) 

where P is a rank-1 projector andr], £ are any pair of real numbers. The eigenvalues 
o/adi are thenztt; (each with multiplicity d—1) andO (with multiplicity d 2 —2d+2). 

Proof. Let Ai > A2 > ■ ■ • > be the eigenvalues of L arranged in decreasing 
order, and let I62), . . . , \bd) be the corresponding eigenvectors. We may assume, 
without loss of generality, that the \b r ) are orthonormal. We have 

adifl&rX&al) = [L, \b r )(b s \] = (A r - A s )|6 r )(6 s | (222) 

So the eigenvalues of ad^ are A r — A s . Since L is not a multiple of the identity we 
must have A r 7^ \ r +i f° r some r in the range 1 < r < d — 1. We then have that 
A s — At 7^ if either s<r<tort<r<s. There are 2r(d — r) such pairs s, t. So 

rank(ad L ) > 2r(d- r) > 2(d- 1) (223) 

Suppose, now that the lower bound is achieved. Then r(d — r) = d — 1, implying 
that r = 1 or d — 1. Also we must have X s = A s+ i for all s/r. So cither 

L = A 2 /+(A 1 -A 2 )|6 1 )(6 1 | (224) 



L = X d -iI-(X d -i-X d )\b d )(b d \ 
Either way L and the spectrum of ad^ are as described. 
The final ingredient needed to complete the proof is 



(225) 
□ 
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Lemma 12. Let L r , C rs t and C r be as defined in the statement of Theorem 7. 
Suppose that the C r have the spectral decomposition 

C r = P r - Pj (226) 

where P r is a rank d — 1 projector which is orthogonal to its own transpose. Let I, 
e' r be as in the statement of Corollary 10. Then there is a fixed sign e = ±1 such 
that 

n r = ee' r L r - (227) 

is a rank-1 projector for all r. 
Proof. Define 

L' r = e' r L r - l —±l (228) 
Then it follows from Corollary 10 that 

Tr(4) = 1 (229) 

for all r, 

Tr(L' r L'J = d ~^^ (230) 

for all r, s, and 

d 2 

]T L' r = dl (231) 

r=l 

It is also easily seen that if we define C' rst = e' r e' s e' t C rst then 

d 2 

[L' r) L' s }=J2C' rst L't (232) 
t=i 



and 

C' r = P' r - P' r T (233) 

where P' r is a rank-1 projector which is orthogonal to its own transpose (see the 
first part of the proof of Theorem 7) . In particular 

rank(ad L; ) = 2(d- 1) (234) 

and the eigenvalues of ad^/ all equal to ±1 or 0. So, taking account of the fact 
that Tr(L' r ) = 1, we can use Lemma 11 to deduce that there is a family of rank-1 
projectors H' r and signs £ r = ±1 such that 

L' r = + (235) 

If C r = +1 (respectively -1) for all r then Eq. (227) holds with n r = U' r and e = +1 
(respectively —1). Also, if d = 2 then L' r is a rank-1 projector irrespective of the 
value of Cr, so Eq. (227) holds with II r = L' r and e = +1. The problem therefore 
reduces to showing that if d > 2 it cannot happen that £ r = +1 for some values 
of r and —1 for others. We will do this by assuming the contrary and deducing a 
contradiction. 

Let m be the number of values of r for which Cr = +1- We are assuming that 
m is in the range 1 < m < d 2 — 1. We may also assume, without loss of generality, 
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that the labelling is such that £ r = +1 for the first m values of r, and —1 for the 
rest. So 

4 = 1? , if? '- m (236) 

Now define 

f rst = Tr(L' r L' s L' t ) (237) 

Eqs. (230) and (231) mean that the same argument which led to Eq. (9) can be 
used to deduce 

L' r L' s = *±i \J2 TrstL'A - K 2 rs I (238) 
Since L[ is a projector it follows that 

L[L' a = (L[) 2 Li s = ^±1 f E^i L ij ^ K ^ ^ 
By essentially the same argument which led to Eq. (118) we can use this to infer 

«> 2 = ITT f > + (5W MM (240) 



where T[ is the matrix with matrix elements T[ rg and ||ei)) is the vector defined by 

2d 
7+1 



Eq. (116). As before | j e i )} is an eigenvector of T[ with eigenvalue ^pj. Consequently 



the matrix 

Qi = ^T[-2|| ei ))« ei || (241) 

is a projector. But that means Tr(Qi) must be an integer. We now use this to 
derive a contradiction. 

It follows from Eq. (236) that 

iK? = {%-. )T d -. T , (242) 



Consequently 



and so 



{K\ r r < m 

2(d-2) d-A ts2 ^ ( 243 ) 

- ^rK( r r>m 



d+l d 



Tr(Q 1 ) = ^^f lrr -2 

= rf+1 _^ + 2^2) 

So if Tr(<5i) is an integer (Ad 2 + 2n(d — 2)) /d 3 must also be an integer. But the 
fact that 1 < m < d 2 , together with the fact that d > 2 means 

4 4g+Md-2) < 2 (245) 
a a- 3 
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If d = 3 or 4 there are no integers in this interval, which gives us a contradiction 
straight away. If, on the other hand, d > 5 there is the possibility 

Ad 2 + 2m(d - 2) 



d 3 

implying 

d 2 (d-A) 



1 (246) 

(247) 



2(d-2) 

This equation has the solution d = 6, m = 9 (this is in fact the only integer solution, 
as can be seen from an analysis of the possible prime factorizations of the numerator 
and denominator on the right hand side). To eliminate this possibility define 

L'r = -,!- L'cP+i-r (248) 



for all r. It is easily verified that 

t 

~ d+1 

d 2 

= dl (250) 



Tr(L^) = (249) 



£ r= a ; „ - A - ( 251 ) 



(252) 



and 

iT r r < d 2 — m 

. \I — n r r > d 2 — m 

So we can go through the same argument as before to deduce 

2 d 2 (d-4) 
a — m = —— — -f 
2(d-2) 

Eqs. (247) and (252) have no joint solutions at all with d ^ 0, integer or otherwise. 

□ 

To complete the proof of Theorem 7 observe that Eqs. (210) and (227) imply 

Tr(II,,n s ) = (253) 

So the II r arc a SIC-set. Moreover 

L r = e r (IT + al) (254) 
where e r = ee' r and a = (el — l)/d. 



6. The Algebra sl(d,C) 

The motivation for this paper is the hope that a Lie algebraic perspective may 
cast some light on the SIC-existence problem, and on the mathematics of SIC- 
POVMs generally. We have focused on gl(d, C) as that is the case where the con- 
nection with Lie algebras seems most straightforward. However, it may be worth 
mentioning that a SIC-POVM also gives rise to an interesting geometrical structure 
in sl(d, C) (the Lie algebra consisting of all trace- zero d x d complex matrices). 



:S2 



Let II r be a SIC-set and define 



So B r g sl(d,C). Let 

(A, A') = Tr(adAadA') = 2dTr(AA') (256) 
be the Killing form [55] on sl(d, C). Then 

(Br,B s ) = \ 1 1 r = S (257) 

t cP-i rz r s 

So the B r form a regular simplex in sl(d, C). Since sl(d, C) is d 2 — 1 dimensional 
the B r are an overcomplete set. However, the fact that 

d 2 

Y B r = (258) 

r=l 

means that for each A g sl(d, C) there is a unique set of numbers a r such that 

d 2 

A = J2a r B r (259) 

r=l 

and 

d 2 

^a r = (260) 

r=l 

The a r can be calculated using 

d 2 - 1 

o r =-^-(A, J B r ) (261) 

Similarly, given any linear transformation M: s\(d, C) — > sl(ti, C), there is a unique 
set of numbers M rs such that 

d 2 

MB r = Y M rs B s (262) 

s=l 

and 

d 2 d 2 

^ M rs = ^ M sr = (263) 

S = l 8=1 

for all r. The M rs can be calculated using 

M rs = —^p—(B s ,MB r ) (264) 

In short, the B r retain many analogous properties of, and can be used in much the 
same way as, a basis. It could be said that they form a simplicial basis. 



33 



7. Further Identities 

In the preceding pages we have seen that there are five different families of ma- 
trices naturally associated with a SIC-POVM: namely, the projectors Q r together 
with the matrices 

Jr = Q r - Qj (265) 
Rr = Q r + Qr (266) 

R r = Q r + Qj + M\e r ))((e r \\ (267) 
T r =^Q r + ^\\e r ))((e r \\ (268) 

(see Section 3). As we noted previously, it is possible to define everything in terms 
of the adjoint representation matrices J r and the rank-1 projectors ||e r ))((e r ||: 

Qr = ~Jr(Jr + I) (269) 
R r = Jr (270) 

R r = J r 2 +4||e r ))((e r || (271) 
Tr = ^YyMJr + + ^yllOXM (272) 

In that sense the structure constants of the Lie algebra, supplemented with the 
vectors ||e r )), determine everything else. 

In the next section we will show that there are some interesting geometrical 
relationships between the hypcrplanes onto which Q r , Qj and R r project. In this 
section, as a preliminary to that investigation, we prove a number of identities 
satisfied by the Q, J and R matries. We start by computing their Hilbert-Schmidt 
inner products: 

Theorem 13. For all r, s 

t In n\ d36 ™ + d 2 -d-l 

Tr{Q r Qs) = ( d+ i)2 ^ 273 ) 



MQrQl) = "-^f (274) 



Tr(J r J s ) = "^f^ (275) 

Tr(RrR s ) = ^-^[^ (276) 

Tr (J r R s ) = (277) 

Proof. Let us first calculate some auxiliary quantities. It follows from the definition 
of T r , and the fact that the matrix P = gG defined by Eq. (63) is a rank d projector, 
that 



Ti(T r T s ) — ^ TruvTs' 
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Also 



^ K uv G ru G us G sv G vr 

u,v=l 

, d 2 ,4 d 2 



- V a 2 a 2 + V p p p p 



d-, 

U=l 'U , U — 1 

d 2 (d<5 rs + d + 2) d 4 , 

(d + 1) 3 ~ + d + i |Pr '' 

d 2 (d£ rs +d + 2) t d 2 k2 



(d+1) 3 d+l 
d 2 (d{d + 2)S rs + 2d + 3) 
(d + l) 3 



Tr (T r Tj) = J2 T ruvT S : 

lt,V = l 

d 2 

^ ^ G ru G su I ^ ^ G uv G uv G vr G v 



2d d 



^ ^ ^ G ru G suG ur Gu 



u=l 

2d 2 



1 + A 2 



(d+l) 2V rs; 
2d 2 (d8 rs + d + 2) 
(d+l) 3 



(278) 



(279) 



where we made two applications of Eq. (23) (i.e. the fact that every SIC-POVM is 
a 2-design). Finally, it is a straightforward consequence of the definitions of T rl Tj 
and || e r )) that 



«e r ||Tje r )) = ((^HTJIK)) 



— V T K 2 K 2 

2^ / j SUV TU TV 

u.v — 1 



1 / 

2d( d I 1) * ^ ^srr H~ ^ ^ ^ T srv -\- d ^ ^ T sur -\- ^ ^ X 1 ^ 



-s // r 

-u— 1 u— 1 u,v— 1 



(3A 2 S + 1) 



2(d + l) 

2(d+l) 2 v ; 
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and 

u— 1 

_ dS rs +d + 2 
= 2(d + l) 
Using these results in the expressions 



(281) 



Tr(Q r Q s ) = Tr( ( ^T r - 2\\e r ))((e r \\] (^T s - 2\\e s ))((e s \\ ) ) (282) 

and 

Tr(Q r Qj) =Trf (^T r -2||e r ))((e r ||) (^-2||* M *ll)l (283) 

the first two statements follow. The remaining statements are immediate conse- 
quences of these and the fact that 

J r = Q r - Qj (284) 
R r = Q r + Qj (285) 

□ 



Now define 

1 d2 

im = -,Eh ( 286 ) 



r=l 

where ||r)) is the basis defined in Eq. (117). The following result shows (among 
other things) that the subspaces onto which the Q r (respectively Qj , R r ) project 
span the orthogonal complement of \\vq)}. 

Theorem 14. For all r 

Q r \\v )) = Qj\\v Q )) = J r \\v )) = R r \\v )) = (287) 



Moreover 



a a j2 

E^ = E^ = ^xt( 7 - K»<M) (288) 



d + i 

r—1 r—1 

d 2 

E Jr = (289) 

r=l 

d2 9d 2 

E^^tt^-inkm) ( 29 °) 

r—1 

Proof. Some of this is a straightforward consequence of the fact that J r is the 
adjoint representative of n r . Since 



En s = dJ (291) 
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we must have 

d 2 d 2 

J rst n t = ad n , n s = (292) 

S,t=l 8=1 

In view of the antisymmetry of the J rst it follows that 

,i 2 



Y Jr = (293) 



and 



Using the relations 



J P ||wo))= (294) 



Q r = \j r {J r +I) (295) 

Qj = \j r {J r -I) (296) 
R r . = J r 2 (297) 



we deduce 



Q r ||«o» = Qj\\vo)) = Rr\\vo)) = (298) 
It remains to prove Eqs. (288) and (290). It follows from Eq. (120) that 

d 2 d 2 d 2 

^Q rst = ^^T rat -2^(( S ||e r ))((e r || t » 

r—1 r—1 r—1 

= {d +l)Kl t - d -±^Y. K rs K2 rt 



d 

r=l 



d 2 S st - 1 
d+l 



(299) 



from which it follows 



r—1 r—1 



Eq. (290) follows from this and the fact that R r — Q r + Qj. 



Geometrical Considerations 



□ 



In this section we show that there are some interesting geometrical relationships 
between the subspaces onto which the operators Q r , Qj and R r project. The 
original motivation for this work was an observation concerning the subspaces onto 
which the R r project. R r is a real matrix, and so it defines a 2(d — 2) subspace 

,2 

in K , which we will denote Ji r . We noticed that for each pair of distinct indices 
r and s the intersection "R r n 3J S is a 1-dimensional line. This led us to speculate 
that a set of hyperplanes parallel to the Jl r might be the edges of an interesting 
polytope. We continue to think that this could be the case. Unfortunately we have 
not been able to prove it. However, it appears to us that the results we obtained 
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while trying to prove it have an interest which is independent of the truth of the 
motivating speculation. 

We will begin with some terminology. Let P be any projector (on either R w 
or C^), let CP be the subspace onto which P projects, and let \ijj) be any non-zero 
vector. Then we define the angle between and CP in the usual way, to be 



(301) 



(so 9 is the smallest angle between \ifj) and any of the vectors in CP). 

Suppose, now, that P' is another projector, and let CP' be the subspace onto 
which P' projects. We will say that CP' is uniformly inclined to CP if every vector in 
CP' makes the same angle 9 with CP. If 6 = this means that CP' C CP, while if 6 = ? 
it means CP' _L CP. Suppose, on the other hand, that < 9 < |. Let \u[), . . . , \u' n ) 
be any orthonormal basis for CP', and define \u r ) = scc9P\u' r ) . Then (u r \u r ) = 1 
for all r. Moreover, if P, P' are complex projectors, 




« + e^u' s \P\u' r + e # <) = 2cos 2 e(l + Re (e** («r|u a ))) 



(302) 



for all <f> and r ^ s. On the other hand it follows from the assumption that CP' is 
uniformly inclined to CP that 



« + e'V.I-PK + e**<) = 2 cos 2 9 



(303) 



for all </> and r ^ s. Consequently 

(u r \u s ) = S rs (304) 

for all r, s. It is easily seen that the same is true if P, P' are real projectors. 

Suppose we now make the further assumption that dim(CP') = dim(CP) = n. Then 
. . . , \u n ) is an orthonormal basis for CP, and we can write 



P = 2J K)(u r | 

r=l 

p' = j2\<)«\ 



Observe that 



(u' r \u s ) — (u' r \P\u s ) — cos9(u r \u s ) = cost 
for all r, s. Consequently 

P'\u r ) = cos9\u r ) 

for all r. It follows that 

ll P 'l^>ll 



cos 6{u r \^)\u' r ) 



= cos ©II |V) 



for all \ip) € CP. So CP is uniformly inclined to CP' at the same angle 
It follows from Eqs. (305) and (306) that 

PP'P = cos 2 9P 

P'PP' = cos 2 9P' 



(305) 
(306) 

(307) 
(308) 

(309) 



(310) 
(311) 



Eq. (310), or equivalently Eq. (311), is not only necessary but also sufficient for 
the subspaces to be uniformly inclined. In fact, let CP, CP' be any two subspaces 
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which have the same dimension n, but which are not assumed at the outset to be 
uniformly inclined, and let P. P' be the corresponding projectors. Suppose 

PP'P = cos 2 OP (312) 

for some 9 in the range < 9 < J . It is immediate that CP = CP' if 9 = 0, and 
CP _L CP' if 6 = ^. Either way, the subspaces are uniformly inclined. Suppose, on 
the other hand, that < 9 < |. Let \ui), . . . , \u' n ) be any orthonormal basis for CP', 
and define \u r ) — sec9P\u' r ) . Eq. (305) then implies 

n n 

P = sec 2 9^P\u' r ){u' r \P = \u r )(u r \ (313) 

r—1 r— 1 

Given any G T we have 

n 

|^)=P|^)=VJ(« r |V;)K) (314) 

r=l 

Since dim(CP) = n it follows that the |u r ) are linearly independent. In particular 

n 

|w r ) = P\u r ) = y~](u a \u r )\u a ) (315) 

s=l 

Since the |w r ) are linearly independent this means 

(u s \u r ) = 5 rs (316) 

So the |u r ) are an orthonormal basis for CP. It follows, that if \ip') is any vector in 
CP', then 

||P|^)|| = 



r=l 



COS 9 



Y^{u' r \4>')\Ur, 



r=l 



= cos $11 IV 7 ') II (317) 



implying that CP' is uniformly inclined to CP at angle 9. 

It will be convenient to summarise all this in the form of a lemma: 

Lemma 15. Let CP, CP' be any two subspaces, real or complex, having the same 
dimension n. Let P , P' be the corresponding projectors. Then the following state- 
ments are equivalent: 

(a) CP' is uniformly inclined to CP at angle 9. 

(b) CP is uniformly inclined to CP' at angle 9. 
(c) 

PP'P = cos 2 9P (318) 

(d) 

P'PP' = cos 2 BP' (319) 
Suppose these conditions are satisfied for some 9 in the range < 9 < ^, and 
let \v,i), . . . \u n ) be any orthonormal basis for CP. Then there exists an orthonormal 
basis \u'i), . . . , \u' n ) for CP' such that 

P'\u r ) = cos0|O (320) 
P\u' r ) = cos 9\u r ) (321) 

We are now in a position to state the main results of this section. Let Q r 
(respectively Q r ) be the subspace onto which Q r (respectively Qj) projects. We 
then have 
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Theorem 16. For each pair of distinct indices r, s the subspaces Q r , Q r have the 
orthogonal decomposition 



whe 



Q° rs J_ Q r 

Q ™ B —I— Q rr 



Q =Q° 

Qr = Q°rs 

dim(Q° s ) = 1 
dim(Q° s ) = 1 



(322) 
(323) 



dim(Q r 
dim(Q r 



= d-2 
= d-2 



We ho 



' l ' 

, d+l, 



1 

/d+l 



(a) Relation of the subspaces Q r and Q s : 

^ 1 ^ g — I — ^ CLThd ^ ^ — I — st 1 " 

(2) Q° s and Q° r are inclined at angle cos 

(3) Q rs and Q sr are uniformly inclined at angle cos -1 

(b) Relation of the subspaces Q r and Q s : 

^ 1 J ^3 g — I — ST* G'fbd S — ^ — ST* ' 

(2) Q° s and Q^ r are inclined at angle cos -1 (jxj) • 

(3) Q rs and Q sr are uniformly inclined at angle cos" 1 ( ^^j j ■ 

(c) Relation of the subspaces Q r and Q s : 

(2) Q° s and Q° r are inclined at angle cos -1 (jrj) ■ 

The relations between these subspaces are, perhaps, easier to assimilate if pre- 
sented pictorially. In the following diagrams the line joining each pair of subspaces 
is labelled with the cosine of the angle between them. In particular a on the line 
joining two subspaces indicates that they are orthogonal. 

On no 




(J 



1 

d+l 




/d+l 



Q° 



Q° 



d 
d+l 
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We will prove this theorem below. Before doing so, however, let us state the other 
main result of this section. Let 5i r be the subspace onto which the R r project. Since 

— ,2 

R r is a real matrix we regard 3? r as a subspace of R . We have 

Theorem 17. For each pair of distinct indices r, s the subspace 5l r has the decom- 
position 

X = Rrs © ffi R rs (324) 
where s , 3£* s} 3? r s are pair-wise orthogonal and 

dim(IR° s ) = 1 dim(0^J = 1 dim(3? rs ) = 2d - 4 (325) 

We have 

(!) S rs = 

(2) ft* s _L % sr and \. s _L $} sr . 

(3) 3£* s and 0Q r are inclined at angle cos _1 (^x)- 



(4) 5i rs and 5l sr are uniformly inclined at angle cos 




In particular, the subspaces 5l r and Jl s intersect in a line. 

In diagrammatic form the relations between these subspaces are 




where, as before, each line is labelled with the cosine of the angle between the two 
subspaces it connects. 

Proof of Theorem 16. Let || 1)) , . . . , ||<i 2 )) be the standard basis for "K^ , as defined 
by Eq. (117). For each pair of distinct indices ?*, s define 

\\f rs ))=iVd + lQ r \\s)) (326) 
II/;,)) = -iVd+lQ?\\a)) (327) 

The significance of these vectors is that ||/ rs ))((/ rs || (respectively ||/ r * s ))((/ r *J) wm 
turn out to be the projector onto the 1-dimensional subspace Q° s (respectively Q° s ). 



41 



Note that the fact that Q r is Hermitian means 

Qj = Q* r (328) 

(where Q* is the matrix whose elements are the complex conjugates of the corre- 
sponding elements of Q r ). Consequently 

= («„»)* (329) 

for all r, s, t. 

It is easily seen that ||/ rs )), ||/* s )) are normalized. In fact, it follows from 
Eqs. (116) and (120) that 

«/r.||/r.» = (d+l)«*||Or||a)) 



d ' 

(d+lf 



2(d+l)« S ||e r »«e r || S » 



(K 2 - K A ) 

^ \ 7*s rs ) 

= 1 (330) 
for all r ^ s. In view of Eq. (329) we then have 

«/Tall/;.» = («/r.ll/r.»)* = l (331) 

for all r s. The fact that Q r Qj = means we also have 

«/„ll/«»= (332) 

for all r =/= s. 

Note that, although we required that r ^ s in the definitions of ||/ rs )), ||/* s )), 
the definitions continue to make sense when r = s. However, the vectors are then 
zero (as can be seen by setting r = s in Eq. (121)). 

The vectors ||/ rs )), ||/* s )) satisfy a number of identities, which it will be conve- 
nient to collect in a lemma: 

Lemma 18. For all r ^ s 

ll/»» = -||/^))+<y|(||e.»-||e r ») (333) 

Il/; S )> = -Il/ S r)}-M/Iflia-Iia) (334) 



d 

(where ||e r )) is the vector defined by Eq. (116)) 

QAU) = \\U) Q?W.)) = (335) 

Orll/r.»=0 QM.))=0 (336) 

QsWU) = -^WU) Q?W.)) = -rf^illO ( 337 ) 

QsWU) = -JTiWfsr)) QsWU) = -^IIO ( 338 ) 

«/rJ/-r» = «/r*sllO = ( 339 ) 
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= ((/; s iio = - ( 34 °) 



Proo/. It follows from Eqs. (116) and (120) that 



((ill/,,)) + «*IIO = «Vd+T(Q rfs - Q srt ) 



d + 1 



= iVcT+T ^ — (T rts - T srt ) 

-2(((t||e r ))((e r || S ))-((r||e s ))((e 8 ||i)) 

= iy|(««-«t||er») (341) 

where we used the fact that T rts = T srt in the third step, and the fact that ((t||e s )) is 
real in the last. This establishes Eq. (333). Eq. (334) is obtained by taking complex 
conjugates on both sides, and using the fact that the vectors ||e s )) are real. 

Eqs. (335) and (336) are immediate consequences of the definitions, and the fact 
that Q r Qj = 0. Turning to the proof of Eqs. (337) and (338), it follows from 
Eqs. (119) and (120) that 

Q.||e a ))=0 (342) 
Using this and the fact that Q s \\f* r )) = in Eq. (333) we find 

QJ/ra» = -*yf<?Jk» (343) 

Since 



\er)) = J^(\\r)) + H))) (344) 



i 2(d+l) 

and taking account of the fact that Q s \\vq)} = (see Eq. (287)) we deduce 



K/^t^W = -^tIIO (345) 

Taking complex conjugates on both sides of this equation we deduce the second 
identity in Eq. (337). 

In the same way, acting on both sides of Eq. (333) with Qj we find 

Qrii/r.»=-io-<Jf«TiM> 



= -IIO-Vl¥T« 



r 



Taking complex conjugates on both sides of this equation we deduce the second 
identity in Eq. (338). 

Turning to the last group of identities we have 

((/-II/-)) = «/r.||Qr||/.r» = ~ ^((U\ U) = (347) 

and 

((U\f: r )) = ((/rJQrlO = -^j((U\U) = (348) 
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The other two identities are obtained by taking complex conjugates on both sides 

of the two just derived. □ 

This lemma provides a substantial part of what we need to prove the theorem. 
The remaining part is provided by 

Lemma 19. For all r ^ s 

QrQsQr = J^jQr ~ J^f\\frsUU\ (349) 

QrQ^Qr = jA^MrMMl (350) 



(d + 1) 

Proof. It follows from Eq. (120) that 

QrQsQr = d -T-QrT s Qr ~ %Qr 1 1 e s )) ((e s 1 1 Q r (351) 

a 

QrQlQr = ^QrTjQr - 2Q r || e s )> «e s || Q r (352) 
In view of Eqs. (344), (287) and the definition of \\f rs )) we have 



Qr\\s))=-i j 1 A irs)) (353) 



Substituting this expression into Eqs. (351) and (352) we obtain 

QrQsQr = ^QrTsQr ~ JJ^W frs)) ((frs\\ (354) 
QrQlQr = —J-QrTjQr ~ II /»» «/r.|| (355) 

The problem therefore reduces to showing 

Q r T s Q r = tttt^Qt (356) 

QrTjQr = jA^\\f r .))ifr.\\ (357) 

Using Eq. (120) we find 

((a\\Q r T s Q r \\b)) = ^^{{a\\T r T s T r \\b}) 
d z 



1 / 2(d + l) 

2 I d 



K? a {(e r \\T s T r \\b}) + K? b {(a\\T r T s \\e r )) 



+ ^±^-KlK 2 rb ((e r \\T s \\e r )) (358) 



((a\\Q r TjQ r \\b)) = ^±fil(( a \\T r TjT r \\b)) 

_ 1 / 2(rf+ 1) 
2\ d 
2{d+l) 



K 2 ra {{e r \\TjT r \\b)) + K 2 M\T r Tj\K))) 
K 2 a K 2 b ((e r \\Tj\\e r )) (359) 
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Using the definitions of T r , \\e r )) and Eq. (23) (the 2-design property) we find, after 
some algebra, 

((a\\T r T s T r \\b)) = j^-^ (K 2 ra T rsb + K 2 rb T ras + K 2 rs T rab + K 2 ra K 2 rb ) (360) 

3 

«e r ||T s T r ||&)) = 2 (^—A 2 (2K 2 s K 2 rb + K 2 rb + T rsb ) (361) 



2( d + 1 ), 

3 

:(a||T r T s ||e r )) = 2 (^yy) ' + Ka + T ras ) (362) 

«e r ||T s ||e r » = {3K 2 S + l) (363) 



and 



({a\\T r TjT r \\b)) = ^—(G ra G as G sb G br 



(d + iy 



K 2 a T rsb + R 2 b T ras + K 2 a K 2 b 



d 2 

(d + l)T ras T rsb 



(d + 1) 



KraT rs b + ^rb^ras + -^ra^rb ) (364) 



( er ||T s T r r ||6)) = 2 ( ) [K 2 s K 2 b + K 2 b + 2T rsb ) (365) 

(a\\T r Tj\\e r )) = 2 ( — A— ) " [K 2 rs K 2 ra + K 2 ra + 2T ras ) 



2{d+l) 

((e r \\Tj\\e r )) = (3K 2 S + l) (367) 

where in deriving Eq. (364) we used the fact that G ra G as G s bGb r = (d+ l)T ras T rsb 

(in view of the fact that r ^ s). Substituting these expressions into Eqs. (358) 

and (359) we deduce Eqs. (356) and (357). □ 

Now define the rank d — 1 projectors 

Qr. = Or-||/r.)>«/r.|| (368) 

Q?. = 0?-H/r*.»«/?all (369) 

and let Q° s , Q rs , Q° s and Q rs be, respectively, the subspaces onto which ||/ rs )) ((/ rs || , 

Qrs' ll/*s))((/*sll an d Qrs project. It is immediate that we have the orthogonal 
decompositions 

Q r = Q° s 8 Q rs (370) 

Q r = Q° rs © Q rs (371) 

Using Lemma 18 we find 

QsrWfrs)) =Qrs\\fsr)) = (372) 
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implying that Q" s _L Q sr and Q rs _L Q" r , and 

\((frs\\f sr ))\ = (373) 

implying that Q° s and Q° r are inclined at angle cos -1 (^j-j-) . Using Lemma 18 
together with Lemma 19 we find 

QrsQ srQrs — QrsQ sQrs 

QrQsQr 1 1 frsf) (ifrs \\ QsQr ~ QrQs \ \ frsj) tf\frs \ \ 

+ «/r.||Qa||/r.»||/r.)>«/r.|| 

= ^- 1 Qr-^- 1 WfrsHU\ 

= -j^—rQrs (374) 

0+1 

which in view of Lemma 15 implies that Q rs and Q sr are uniformly inclined at angle 
cos -1 (^==). This proves part (a) of the theorem. Parts (b) and (c) are proved 
similarly. 



Proof of Theorem 17. Define 

IISra» = ^(ll/;.)> + ll/r.») (375) 
llflra» = ^(IIO-H/r.») (376) 

By construction the components of ||ft s )), lifts)) m the standard basis are real, so 
we can regard them as <S R d . They are orthonormal: 

((ftsllfts)) = ((ftsllfts)) = 1 and ((ftsllfts)) = (377) 
It is also readily verified, using Lemma 18, that 

Rrhrs)) = \\9rs)) (378) 

Rr\\9rs)) = ||5r.)) (379) 

So 

R rs =Rr- lifts)) ((fts|| ~ ||fts)) ((fts || (380) 

is a rank 2d — 4 projector. If we define 3?° s , 3^ s and 3l rs to be, respectively, 
the subspaces onto which lifts)) ((fts II j lifts)) ((.fts I an( l R rs project we have the 
orthogonal decomposition 

X r = © $i. s © K rs (381) 
It follows from Eqs. (333) and (334) that 

lifts)) = -lift,)) (382) 

implying that 3?° s = 3?° r for all r ^ s. It is also easily verified, using Lemma 18, 
that 

|((fta||ftr))| = J^ (383) 

from which it follows that 3^ s and $} sr are inclined at angle cos~ 1 (^j). We next 
observe that 

R rs = Q rs + Ql (384) 
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Using Lemma 18 once again we deduce 

RrsUsr)) =Rsr\\9rs)) = (385) 

from which it follows that _L 3l sr and Ji rs _L 3?J. r . Finally, we know from 
Theorem 16 that Qj s Q sr = Q rs Qj r — 0. Consequently 



R-rsRsrR-rs QrsQ srQrs QrsQ srQ 

d n i d ^ T 



rs 



d + i^ rs d + i^ rs 

In view of Lemma 15 it follows that Ji rs and Jl sr are uniformly inclined at angle 

Further Identities. We conclude this section with another set of identities in- 
volving the vectors ||/ rs )), ||/ r * s )), \\g rs )) and \\g rs }}. 
Define 

2d „ d + 1 ,, .. . „ 

where \\vo)} is the vector defined by Eq. (286). It is readily verified that 

((e r ||e r ))=0 and ((e r \\v )) = (388) 

So ||e r )), \\vq)) is an orthonormal basis for the 2-dimensional subspace spanned by 
|| e r )), \\v )}. Note that 

Q r ||e r )) = Qj\\e r )) = R r \\e r )) = (389) 

We then have 
Theorem 20. For all r 

d 2 

d+~ 



E II /-))((/- II =Qr (390) 
s=l 

* 

E \\frs)){(f: s \\=Qr (391) 

8=1 
r, d 2 

J^l E \\9r S ))((9r S \\=Rr (392) 



s=l 



2 d 



I E \\9r.)){{3r,\\=Rr (393) 



s—1 



and 

,1 



~Y E \\f«-UU\=Qr +\\e r ner\\ + ¥ ^ Y (l-\\v Uv \\) (394) 



d 

s=l 
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2 



<p 

E 


IIO((/a*r 


s=l 
( s /r) 




d 2 




E 




s=l 




d 2 




E 


\\9sr)){(gsr 


s=l 
( a #r) 





(395) 



= i? r (396) 
5 , 4(d-l) M .... „ , 4 



J-IN»««o| 

(397) 



Proof. It follows from the definition of ||/ rs )) that 

d 2 ci 2 
^7 E ll/«»«/r.|| = E Qr\\4((4Qr 

8=1 8=1 

= Q r (e 11*11 j 2r 

= Q r (398) 

where in the second step we used the fact that Q r ||r)) = (as can be seen by setting 
r = s in Eq. (121)). Eq. (391) is obtained by taking the complex conjugate on both 
sides. 

We also have 

d 2 d 2 
J-7 E H/r.»«#.ll = - E Orll*»«*IIOr 

8=1 8=1 

= ~Qr [Ell«»«4 ] 0? 



s=l 
QrQr 

= (399) 



Taking the complex conjugate on both sides we find 



-t^t E H#.»«/rJ= (400) 

(s^r) 

Consequently 

d 2 d 2 
r E \\9r.U9r.\\ = (ll/«»«/r.|| + II /r*.» ((frs II 



S=l S— 1 



+ ll/r.»«^.ll + ll/?.»«/r, 

i? r (401) 
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Eq. (393) is proved similarly. 

To prove the second group of identities we have to work a little harder. Using 
Eqs. (116) and (120) we find 



E H/-»«/-H & » = ^EH^Ik»((r||Q s ||&}) 



d- - 

s=l s = l 



d 2 (d- 1) Z^V ar fc sa sr ; 



srb 



- Ki r Ki b T sar + Ki a K%.Ki b ) (402) 

(where we used the fact that <5 s ||s)) = in the first step). After some algebra we 
find 



TrlJ (403) 



T r6a ) (404) 



d+1 



T rba ) (405) 



d+1 



\/d 








d- 


1-2 



+ (d^ d - + (dTTFj (406) 

where we used Eq. (23) to derive the first expression. Substituting these expressions 
into Eq. (402) and using 



(407) 

we deduce Eq. (394). Taking complex conjugates on both sides we obtain Eq. (395). 

Eq. (396) is an immediate consequence of Eq. (392) and the fact that \\g sr )) = 
-\\g rs )) for all r,s. 
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To prove Eq. (397) observe that it follows from Eqs. (394)-(396) that 

d 2 d 2 

E (iiora + iio«/ji) = E ( 2 ib-»«<u - iio«/ji - uora 



S=l 8=1 



(^.-(d-lJUerK^H 

-^(l-\\vo)){{vo\\)) (408) 



Hence 



E iifl.r)>«^ii = E + no(a 



d-3 ^ rf_3 

s— 1 s— 1 



-iiora - no 



^ + ^r"^^" + ( d + i) 4 ( d -3) (- 



(409) 
□ 



9. The P-P t Property 

In the preceding sections the Q-Q T property has played a prominent role. In 
this section we show that in the particular case of a Weyl-Heisenberg covariant SIC- 
POVM, and with the appropriate choice of gauge, the Gram projector (defined in 
Eq. (63)) has an analogous property, which we call the P-P T property. Specifically 
one has 

PP T = P T P = \\h)) ((h\\ (410) 

where \\h}} is a normalized vector whose components in the standard basis are all 
real. In odd dimensions the components of \\h)) in the standard basis can be simply 
expressed in terms of the Wigner function of the fiducial vector. It could be said 
that the projectors P and P T are almost orthogonal (by contrast with the projectors 
Q r and Qj which are completely orthogonal). More precisely P has the spectral 
decomposition 

P = P+\\h))((h\\ (411) 
where P is a rank (d — 1) projector with the property 

PP T = (412) 

This means that the matrix 

J P = P - P T (413) 

is a pure imaginary Hermitian matrix with the property that Jj, is a real rank 
2d — 2 projector (c.f. the discussion in Section 4). 

Although we are mainly interested in the P-P T property as it applies to SIC- 
POVMs, it should be noted that it actually holds for any Weyl-Heisenberg covariant 
POVM (with the appropriate choice of gauge) . So we will prove the above propo- 
sitions for this more general case. 



Let us begin by fixing notation. Let |0), . . . , \d — 1) be an orthonormal basis for 
d-dimcnsional Hilbert space and let X and Z be the operators whose action on the 
\r) is 

X|o) = |o + l) (414) 
Z\a)=U) a \d) (415) 

where ui = e^r and the addition of indices in the first equation is mod d. We then 
define the Weyl-Heisenberg displacement operators by (adopting the convention 
used in, for example, ref. [16]) 

D p = t PiP2 X Pi Z P2 (416) 

where p is the vector (pi,P2) (pi, P2 being integers) and r = e^d 1 . Generally 
speaking the decision to insert the phase t PiP2 is a matter of convention, and many 
authors define it differently, or else omit altogether. However, for the purposes of 
this section it is essential, as a different choice of phase at this stage would lead to a 
different gauge in the class of POVMs to be defined below, and the Gram projector 
would then typically not have the P-P T property. 

Note that t 2 = r d — oj in every dimension. If the dimension is odd we can write 
t = . So r is a d th root of unity. However, if the dimension is even r d = — 1. 
This has the consequence that 

A>+rfu = (-ir p2+u2pi £p (417) 

So in even dimension p = q (mod d) does not necessarily imply D p = _D q (although 
the operators are, of course, equal if p = q (mod 2d)) 
In every dimension (even or odd) we have 

Dl = D_ p (418) 

for all p 

(D p ) n = D np (419) 

for all p, n and 

D p D (1 = t^D p+ci (420) 
for all p, q. In the last expression (p, q) is the symplectic form 

(P,q) =P2?i -Pi<?2 (421) 

Now let \tp) be any normalized vector (not necessarily a SIC-fiducial vector), and 
define 

\%)=D v \ip) (422) 

Let 

It is easily seen that 
for all p. 



£=Ei^>^pi ( 423 ) 



pG2 



[D P7 L]=Q (424) 
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We now appeal to the fact that there is no non-trivial subspace of Jf^ which 
the displacement operators leave invariant. To see this assume the contrary. Then 
there would exist non-zero vectors \<f>), \x) such that 

(<P\D p \ X ) = (425) 

for all p. Writing the left-hand side out in full this gives 

d-l 

]r^ Q (0|a + Pl )Hx)=O (426) 

a=0 

for all pi,P2- Taking the discrete Fourier transform with respect to pi , we have 

(<t>\a + Pl )(a\ X ) = (427) 

for all a, pi, implying that cither |</>) = or \x) = — contrary to assumption. We 
can therefore use Schur's lemma [55] to deduce that 

L = kl (428) 

for some constant k. Taking the trace on both sides of this equation we infer 
that k — d. We conclude that is a POVM. We refer to POVMs of this 

general class as Wcyl-Hciscnbcrg covariant POVMs. We refer to the vector \ip) 
which generates the POVM as the fiducial vector (with no implication that it is 
necessarily a SIC-fiducial). 

Now consider the Gram projector 

P= E ^qllP»«q|l (429) 

p,q G Z5 

where 

P P,q = ^ P IV>q) (430) 

and where we label the matrix elements of P and the standard basis kets with the 
vectors p, q rather than with the single integer indices r, s as in the rest of this 
paper. We know from Theorem 1 that P is a rank d projector. 
In view of Eqs. (418) and (420) we have 



(p||P||q)) = P p , q 

i r -(p.i)(^|D q _ p |^) 



d 

-l 



a=0 



Hence 

((p||PP T ||q)> = V((p||P||u»((q||P||u); 



u6Z d 



d-l 

r PlP2+9l92, ,U2(u 1 +a+b)-(u 1 +a)p2-(u 1 +b}q 2 

d 2 

a,b,ui ,U2- 



x (ip\a + u 1 + -qi)(a\ip)(b\ilj} 
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d-1 

— _ r PlP2+<?ig2 w P2b+92a^| 
^ a,b=0 

= ((p\\h)){(h\\q)) 
where \\h}} is the vector with components 

d-l 

(( P \\h)) = -=Y, TP1P2ujP2a ^~ a -p^( a ^ ( 433 ) 

v a— 

It is easily verified that \\h)) is normalized, and that ((p\\h)} is real. 

Finally, suppose that the dimension is odd. Then the Wigner function of the 
state is [56, 57] 

W(p) = ~{i>\D p U P Dl\i>) = ^\D 2p U P \i>) (434) 

where Up is the parity operator, whose action on the standard basis is Up \a) = | —a) . 
It is straightforward to show 

((p\\h)) = ^/dW{-2- l p) (435) 

where 2 _1 = (d+ l)/2 is the multiplicative inverse of 2 considered as an element of 
Zd'. i.e. the unique integer < m < d such that 2m = 1 (mod d). 

10. Conclusion 

A curious fact about SIC-POVMs is that, although they are characterized by 
their being highly symmetric structures, they do not wear this property on their 
sleeve (so to speak). If one casually inspects the components of a SIC-fiducial, 
without knowing in advance that that is what they are, there does not seem to be 
anything special about them at all. Indeed, so far from there being any obvious 
pattern to the components, they seem, to a casual inspection, like a completely 
random collection of numbers. Moreover, this is just as true of an exact fiducial as 
it is of a numerical one (see, for instance, the tabulations in Scott and Grassl [46]). 
It is only when one looks at them through the right pair of spectacles, and takes the 
trouble to calculate the overlaps Tr(II r n s ), that the symmetry becomes apparent. 
The situation is a little reminiscent of a hologram, which only takes on the aspect of 
a meaningful image when it is viewed in the right way. If one wanted to summarize 
the content of this paper in a nutshell it could be said that we have exhibited some 
other pairs of spectacles — other ways of looking at a SIC — which cause its inner 
secrets (or at any rate some of its inner secrets) to become manifest. 

Rather than focusing on the SIC-vectors \if) r ), as is usually done, we have focused 
on the angle tensors 8 rs and 9 rst , and on the T, J and R matrices defined in 
terms of them. This is an important change of emphasis because, rather than 
being tied to any particular SIC, these quantities characterize entire families of 
unitarily equivalent SICs. Like the components of a SIC-fiducial, the angle tensors 
appear, to a casual inspection, like a random collection of numbers. However, if 
one examines the spectra of the T, J and R matrices one realizes that, underlying 
the appearance of randomness, there is a high degree of order. If one then goes on 
to examine the geometrical relationships between the subspaces onto which the Q, 
Q T and R matrices project, as we did in Section 8, one finds yet more instances 
of structure and order. To our minds what is particularly interesting about all of 



-b- Pl ){b\Tp)^\-a-qi){a\i>) 

(432) 
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this is that none of it is obviously suggested by the defining property of a SIC, that 
Tr(II r n 8 ) = l/(d + 1) for r ^ s. 

In the course of this paper we have several times expressed the hope that the 
Lie algebraic perspective on a SIC will lead to a solution to the existence problem. 
Of course, that is only a hope, and it may not be fulfilled. However, we feel on 
rather safer ground when we suggest that the solution is likely to come, if not from 
this investigation, then from one which is like it to the extent that it focuses on a 
feature of a SIC which is not immediately apparent to the eye. 

Specializing to the case of a Weyl-Heisenberg covariant SIC, a fiducial vector 
is a solution to the equations 

\(*\DM* = ^L±± (436) 

Allowing for the arbitrariness of the overall phase of and taking \ip) to be nor- 
malized, this gives us d 2 — 1 conditions on only 2d — 2 independent real parameters. 
The equations are thus over-determined, and very highly over-determined when d 
is large. Nevertheless, they have turned out to be soluble in every case which has 
been investigated to date. It seems likely that progress will depend on finding the 
structural feature which is responsible for this remarkable fact. The motivation for 
this paper is the belief that it may be structural features of the Lie algebra gl(d, C) 
which are responsible. That suggestion may or may not be correct. But if it turns 
out to be incorrect, the amount of effort which has been expended on this problem 
over a period of more than ten years, so far without fruit, suggests to us that the 
solution will depend on finding some other structural feature of a SIC, which is not 
obvious, and which has hitherto escaped attention. 
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